home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / sgesv.z / sgesv
Encoding:
Text File  |  2002-10-03  |  3.9 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGEEEESSSSVVVV((((3333SSSS))))                                                            SSSSGGGGEEEESSSSVVVV((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGESV - compute the solution to a real system of linear equations A * X =
  10.      B,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SGESV( N, NRHS, A, LDA, IPIV, B, LDB, INFO )
  14.  
  15.          INTEGER       INFO, LDA, LDB, N, NRHS
  16.  
  17.          INTEGER       IPIV( * )
  18.  
  19.          REAL          A( LDA, * ), B( LDB, * )
  20.  
  21. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  22.      These routines are part of the SCSL Scientific Library and can be loaded
  23.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  24.      directs the linker to use the multi-processor version of the library.
  25.  
  26.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  27.      4 bytes (32 bits). Another version of SCSL is available in which integers
  28.      are 8 bytes (64 bits).  This version allows the user access to larger
  29.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  30.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  31.      only one of the two versions; 4-byte integer and 8-byte integer library
  32.      calls cannot be mixed.
  33.  
  34. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  35.      SGESV computes the solution to a real system of linear equations A * X =
  36.      B, where A is an N-by-N matrix and X and B are N-by-NRHS matrices.
  37.  
  38.      The LU decomposition with partial pivoting and row interchanges is used
  39.      to factor A as
  40.         A = P * L * U,
  41.      where P is a permutation matrix, L is unit lower triangular, and U is
  42.      upper triangular.  The factored form of A is then used to solve the
  43.      system of equations A * X = B.
  44.  
  45.  
  46. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  47.      N       (input) INTEGER
  48.              The number of linear equations, i.e., the order of the matrix A.
  49.              N >= 0.
  50.  
  51.      NRHS    (input) INTEGER
  52.              The number of right hand sides, i.e., the number of columns of
  53.              the matrix B.  NRHS >= 0.
  54.  
  55.      A       (input/output) REAL array, dimension (LDA,N)
  56.              On entry, the N-by-N coefficient matrix A.  On exit, the factors
  57.              L and U from the factorization A = P*L*U; the unit diagonal
  58.              elements of L are not stored.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSGGGGEEEESSSSVVVV((((3333SSSS))))                                                            SSSSGGGGEEEESSSSVVVV((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      LDA     (input) INTEGER
  75.              The leading dimension of the array A.  LDA >= max(1,N).
  76.  
  77.      IPIV    (output) INTEGER array, dimension (N)
  78.              The pivot indices that define the permutation matrix P; row i of
  79.              the matrix was interchanged with row IPIV(i).
  80.  
  81.      B       (input/output) REAL array, dimension (LDB,NRHS)
  82.              On entry, the N-by-NRHS matrix of right hand side matrix B.  On
  83.              exit, if INFO = 0, the N-by-NRHS solution matrix X.
  84.  
  85.      LDB     (input) INTEGER
  86.              The leading dimension of the array B.  LDB >= max(1,N).
  87.  
  88.      INFO    (output) INTEGER
  89.              = 0:  successful exit
  90.              < 0:  if INFO = -i, the i-th argument had an illegal value
  91.              > 0:  if INFO = i, U(i,i) is exactly zero.  The factorization has
  92.              been completed, but the factor U is exactly singular, so the
  93.              solution could not be computed.
  94.  
  95. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  96.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  97.  
  98.      This man page is available only online.
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.